ON THE SPLITTING METHOD FOR SCHRODINGER-LIKE 
EVOLUTION EQUATIONS 



ZDZISLAW BRZEZNIAK AND ANNIE MILLET 

Abstract. Using the approach of the sphtting method developed by I. Gyongy 
and N. Krylov for parabohc quasi hnear equations, we study the speed of conver- 
gence for general complex-valued stochastic evolution equations. The approxima- 
tion is given in general Sobolev spaces and the model considered here contains both 
the parabolic quasi-linear equations under some (non strict) stochastic parabolic- 
ity condition as well as linear Schrodinger equations 



1. Introduction 

Once the well-posedeness of a stochastic differential equation is proved, an im- 
portant issue is to provide an efficient way to approximate the unique solution. The 
aim of this paper is to propose a fast converging scheme which gives a simulation 
of the trajectories of the solution on a discrete time grid, and in terms of some spa- 
tial approximation. The first results in this direction were obtained for Stochastic 
Differential Equations and it is well-known that the limit is sensitive to the ap- 
proximation. For example, the Stratonovich integral is the limit of Riemann sums 
with the mid-pint approximation and the Wong Zakai approximation also leads to 
Stratonovich stochastic integrals, and not to the Ito ones in this finite dimensional 
framework. There is a huge literature on this topic for stochastic PDEs, mainly 
extending classical deterministic PDF methods to the stochastic framework. Most 
of the papers deal with parabolic PDFs and take advantage of the smoothing ef- 
fect of the second order operator; see e.g. [9], [21], [H], [12], [13], [13], [18] and the 
references therein. The methods used in these papers are explicit, implicit or Crank- 
Nicholson time approximations and the space discretization is made in terms of finite 
differences, finite elements or wavelets. The corresponding speeds of convergence are 
the "strong" ones, that is uniform in time on some bounded interval [0,T] and with 
various functional norms for the space variable. Some papers also study numerical 
schemes in other "hyperbolic" situations, such as the KDV or Schrodinger equations 
as in [3], [3] and [1]. Let us also mention the weak speed of convergence, that is 
of an approximation of the expected value of a functional of the solution by the 
similar one for the scheme obtained by [1] and [6]. These references extend to the 
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infinite-dimensional setting a very crucial problem for finite-dimensional diffusion 
processes. 

Another popular approach in the deterministic setting, based on semi-groups the- 
ory, is the splitting method which solves successively several evolution equations. 
This technique has been used in a stochastic case in a series of papers by I. Gyongy 
and N. Krylov. Let us especially mention reference [10] which uses tools from [20] . 
[T5] and [16], and provides a very elegant approach to study quasilinear evolution 
equations under (non strict) stochastic parabolicity conditions. In their framework, 
the smoothing effect of the second operator is exactly balanced by the quadratic 
variation of the stochastic integrals, which implies that there is no increase of space 
regularity with respect to that of the initial condition. Depending on the num- 
ber of steps of the splitting, the speed of convergence is at least twice that of the 
classical finite differences or finite elements methods. A series of papers has been 
using the splitting technique in the linear and non-linear cases for the deterministic 
Schrodinger equation; see e.g. [1], [Z], [IS] and the references therein. 

The stochastic Schrodinger equation studies complex-valued processes where the 
second order operator iA does not improve (nor decay) the space regularity of the 
solution with respect to that of the initial condition. Well-posedeness of this equation 
has been proven in a non-linear setting by A. de Bouard and A. Debussche [1] ; these 
authors have also studied finite elements discretization schemes for the corresponding 
solution under conditions stronger than that in [3]. 

The aim of this paper is to transpose the approach from [10] to general quasilinear 
complex- valued equations including both the "classical degenerate " parabolic setting 
as well as the quasilinear Schodinger equation. Indeed, the method used in [10] 
consists in replacing the usual splitting via semi-groups arguments by the study of 
p-th moments of — Z'^ where Z^ and Z^ are solutions of two stochastic evolution 
equations with the same driving noise and different families of increasing processes 
Vq and V{ for r = 0, 1, ■ ■ ■ ,di driving the drift term. It does not extend easily to 
nonlinear drift terms because it is based on some linear interpolation between the 
two cases Vq and V{. Instead of getting an upper estimate of the p-th moment of 
Z^ — Z^ in terms of the total variation of the measures defined by the differences 
Vq — V{, using integration by parts they obtain an upper estimate in terms of the 
sup norm of the process {¥^(1) — V{(t),t G [0,T]). 

We extend this model as follows: given second order linear differential operators 

, r = 0, ■ ■ ■ ,di with complex coefficients, a finite number of sequences of first 
order linear operators S\ I > 1 with complex coefficients, a sequence of real- valued 
martingales M', / > 1 and a finite number of families of real- valued increasing 
processes V[, i G {0, 1}, r = 0, ■ ■ ■ , di, we consider the following system of stochastic 
evolution equations 



dZi{t) = Lr{t, ■)Zi{t)dV[{t) + J2 Si{t, ■)Z,{t)dM\t), z = l,--- ,d, 

r I 
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with an initial condition Zi{0) belonging to the Sobolev space if™'^ for a certain 
m > 0. Then under proper assumptions on the various coefficients and processes, 
under which a stochastic parabolicity condition (see Assumptions (Al)- (A4(m,p)) 
in section [2]), we prove that for p G [2, oo), we have 

e( sup \\Z^{t)-Zom'rn) <c(E\\Z,{0)-Zo{OWrn + AA, (1.1) 

where A = sup^ sup^gjQj.] max^ \ V{{t) — VQ{t)\. When the operator Lr = iA + U 
for certain ffist order differential operator L^, we obtain the quasilinear Schrodinger 
equation. Note that in this case, the diffusion operators are linear and cannot 
contain ffist order derivatives. 

As in [To], this abstract result yields the speed of convergence of the following 
splitting method. Let r„ = {iT/n,i = 0, ■ ■ ■ ,n} denote a time grid on [0,T] with 
constant mesh 6 = T/n and define the increasing processes At{n) and Bt{n) = 
At+s{n), where 

. , / M for t G [2k6, {2k + 1)6], 

~ \t-{k + l)6 for t G [(2A; + 1)6, {2k + 2)6]. 

Given a time-independent second order differential operator L, ffist order time- 
independent operators 5' and a sequence {W'', / > 1) of independent one- dimensional 
Brownian motions, let Z, Zn and Cn be solutions to the evolution equations 

dZ{t) = LZ{t)dt + ^SiZ{t)o dWl, 
dZr,{t) = LZ^{t)dAt{n) + SiZn{t) o 

dCn{t) = LCn{t)dBt{n) + Sl{t, ■)Cn{t) O 

I 

where odWj: denotes the Stratonovich integral. The Stratonovich integral is known 
to be the right one to ensure stochastic parabolicity when the differential operators 
Si contain ffist order partial derivatives (see e.g. [ID])- Then (n{2k6,x) = Z{k6,x), 
while the values of Zn{2k6,x) are those of the process Z„ obtained by the following 
spitting method: one solves successively the correction and prediction equations on 
each time interval [iT/n, {i + l)T/n)]: dvt = Lvtdt and then dvt = X]; S^t, ■)vtodWj:. 

Then, one has A = CT/n, and we deduce that E^supjg[o;r] ll^(^) ~ ^n{t)\\^'^ < 

Cn~P. As in [TU], a fc-step splitting would yield a rate of convergence of Cn^''^. 

The paper is organized as follows. Section [2] states the model, describes the evo- 
lution equation, proves well-posedeness as well as apriori estimates. In section [3] 
we prove (11. ip first in the case of time-independent coefficients of the differential 
operators, then in the general case under more regularity conditions. As explained 
above, in section H] we deduce the rate of convergence of the splitting method for 
evolution equations generalizing the quasi-linear Schrodinger equation. The speed 
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of convergence of the non-linear Schrodinger equation will be addressed in a forth- 
coming paper. As usual, unless specified otherwise, we will denote by C a constant 
which may change from one line to the next. 

2. Well-posedeness and first apriori estimates 

2.1. Well-posedeness results. Fix T > 0, F = (J't, t G [0, T]) be a filtration on the 
probability space {Q, J-", P) and consider the following C-valued evolution equation 
on the process Z(t, x) = X(t, x) + iY(t, x) defined for t G [0, T] and x G M.'^: 

di 

dZ{t, x) = Y^ [LrZ{t, x) + Fr{t, x)] dV[ + ^ [SiZ{t, x) + Gi{t, x)] dM^, (2.1) 

Z(0,x) = Zo{x) = Xoix)+tYo{x), (2.2) 

where di is a positive integer, t G [0, T]), r = 0, 1, ■ ■ ■ ,di are real-valued increas- 
mg processes, {Ml,t G [0,T]), / > 1, are independent real-valued {J^t,t G [0,T])- 
martingales, Lr (resp. Si) are second (resp. first) order differential operators defined 
as follows: 



d d 

LrZ{t,x) = Dk(^[ai-\t,x) + ib>,'\t)]D,Z{t,x)'^ + ^ a^,(t, x)D,Z(t, a;), 
j,k=i j=i 

+ [ttrit, x) + ibr{t, x)~\ Z(t, x), (2.3) 
d 

SiZ{t,x) = Y,al{t,x)D,Z{t, x) + [(7i{t,x) + iTi{t,x)\Z{t,x). (2.4) 

Let m > be an integer. Given C-valued functions Z{.) = X{.) + iY{.) and 
((.) = ^(.) + iry(.) which belong to L'^iR'^), let 



(Z,C) = (^,C)o:= / Re{Z{x)C{x))dx = [X{x)ax) + Y{xUx)]dx. 

Thus, we have (X,^) = J^^ X{x) ^{x) dx, so that (Z, C) = (X,0 + (Y,7]). For 
any multi- index a = (ai, 02, ■ ■ ■ , dd) with non- negative integer components ctj, set 
|a| = X]j '^j ^ regular enough function / : M*^ — M, let -0°"/ denote the 

partial derivative (gf^)"^ ■ ■ ■ (af^)"'*(/)- For ^ = 1, ■ ■ ■ , c^, let Dfc/ denote the partial 
derivative For a C-valued function F = Fi + iF2 defined on W^, let D^F = 
D'^Fi + iD°'F2 and -D^-F = -D/c-f^i + iF)kF2. Finally, given a positive integer m, say 
that F G if™ if and only if Fi and F2 belong to the (usual) real Sobolev space 
if™ = if*"'^. Finally, given Z = X + lY and ( = ^ + irj which belong to set 

(Z,C).„= Yl I Re{D'^Z{x)D%{x))dx (2.5) 

a:U<|Q|<m 
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= J2 f [D'^X{x)D''^{x) + D''Y{x)D''r]{x)]dx, 

a:0<\a\<m 

\\Z\\l = {Z,Z)m= V] [ Re{D''Z{x)D-Z{x))dx. (2.6) 

a:\a\<m 

We suppose that the following assumptions are satisfied: 

Assumption (Al) For r = 0, ■ ■ ■ ,di, {V^^,t G [0,T]) are predictable increasing 
processes. There exist a positive constant K and an increasing predictable process 
{Vt,te [0,T]) such that: 

Vq = V^' = 0, r = 0, ■ ■ ■ du Vt<K a.s., (2.7) 

di 

'y^^dV^ + d{M^)t < dVt a.s. in the sense of measures. (2-8) 

r=0 1>1 

Assumption (A2) 

(i) For r = 0,1,- --di, the matrices {ai:''(t,x),j,k = l,---d) and {hi:^{t),j,k = 
1, - ■ ■ d) are (J-'j) -predictable real- valued symmetric for almost every uj,t& [O5 ^] ^iid 
X e M"*. 

(ii) For every t G (0, T] and x, ?/ G M'^ 

J2y'y' [2 E ^) ^^t' - E ^)^'(^' > (2-9) 

j,k=l r=0 1>1 

a.s. in the sense of measures. 

Assumption (A3(m)) There exists a constant K{m) such that for all j, k = 1^ - ■ ■ d, 
r = 0, ■ ■ ■ , (ii, / > 1, any multi-indices a (resp. /3) of length |a| < m + 1 (resp. 
< m), and for every (t,x) G (0,T] x R"^ one has a.s. 

\D''al'^{t,x)\ + |6^''=(t)| + \D''al{t,x)\ + |Z^V(t,a;)| + |D^6,(t,x)| < k{m), (2.10) 

|D°a/(t,x)| + \D''(Tj{t,x)\ + |D"rj(t,x)| < K(m). (2.11) 

Assumption (A4(m,p)) Let p G [2, +00); for any r = O,--- / > 1, the 

processes Fr{t,x) = Fr^i(t,x) + iFr^2(t,x) and Gi(t,x) = Gi^i(t,x) + iGi^2{t,x) are 
predictable, Fr{t, ■) G iy"* and G,.(t, ■) G Furthermore, if we denote 



^m{t) = I [El|i^.(s)|lLrfK'' + Ell^'(^)ll-+l^(^')^]' (2.12) 

then 

E(||Zo||^ + A'l(T)) < +CX). (2.13) 

The following defines what is considered to be a (probabilistically strong) weak 
solution of the evolution equations (12.ip - (]2.2p . 
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Definition 1. A C-valued {J^t) -predictable process Z is a solution to the evolution 
equation fl2.ip with initial condition Zq if 



F / \\Z(s)\\Us < +00] =1, E \Z(sWdK< 00 



>(^^ \\Z{s)\\lds < +00 









E [ 1 




Jo 



dVJs) 



and for every t G [0, T] and $ = + itp, where cf) and tp are C°° functions with 
compact support from to M, one has a.s. 

(Z(t),$) = (Z(0),$) + f; A - E {[ans,■) + ^K•\s)]D,Z{s,■),D,<l>] 

r=0 "'O j,k=l 

d 

+ ■)D,Z{S, ■) + K(S, ■) + lhr{s, ■)]Z{S, ■) + ■), ^] 

+ 5^ / ■)) + Gi{s, ■), $) dMl (2.14) 

Theorem 1. Let m > 1 be an integer and suppose that Assumptions (A1),(A2), 
(A3(m)) and (A4(m,2)) (i.e., forp = 2) are satisfied. 

(i) Then equations (12. ip and (12.21) have a unique solution Z, such that 

e( sup \\Z{t,-)\\l) <Ce(\\Z,\\1 + K^{T)) <oo, (2.15) 
^te[o,T] ^ ^ ^ 

for a constant C that only depends on the constants which appear in the above 
listed conditions. Almost surely, Z e C([0, T], if™"^) and almost surely the map 
[0,T] 3 t^ Z{t, ■) G H"^ is weakly continuous. 

(a) Suppose furthermore that Assumption (A4(m,p)) holds for p G (2, +00). 
Then there exists a constant Cp as above such that 

e( sup \\Z{tr)\\l) <cM\\Z4l + Ki:('iT)). (2-16) 

Proof Set n = i/™, V = H'^+^ and V = H'^-K Then V C C V is a Gehand 
triple for the equivalent norm |(/ — A)™'/^n|i2 on the space Hm- Given Z = X + iY G 
V and C = ^ + ir] eV set 

{Z, Om = {X, Om + {Y, v)m, and (Z, C) := {Z, Qo, 

where {X,^)m and {Y,ri)m denote the duality between the (real) spaces if^+i and 
ffrn-i Pqj, gvery multi-index a, let 

X(a) = {(/3,7) : « = /3 + 7, I7I e {0, ■ ■ ■ , |a|} }. 

To ease notations, we skip the time parameter when writing the coefficients ar, br, 
ai and r/. Then for / > 1, using the Assumption (A3(m)), we deduce that if 
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Z = X + iY ^ iJ™"*"^, we have for |a| < m, 

d 

D"[Si{Z)] = Y,cri{x)D,D"Z+ CiiP,l)D^Z, (2.17) 

with functions Ci{f3, 7) from R'^ to C such that sup;>]^ sup^^gi^d \Ci{f3, 7)(x)| < +00. A 
similar computation proves that for every multi-index a with \a\ < m, r = 0, ■ ■ ■ ,di 

d 

D^[Lr{Z)] = LriD-Z) + ^ Yl Dk{D^al:''D,D^z) + 7)^''^, 

j,fc=l(/3,7)GX(a):|7|=l |/3|<m 

(2.18) 

for some bounded functions Cr{P, 7) from M.'^ into C. Hence for every r = 0, ■ ■ ■ , di, 
one has a.s. : V x f2 — )■ V and similarly, for every I > 1, a.s. Si : V x Q ^ Ti. 
For every A > and Z = X + iY & Ti, let us set 

UaZ := L^Z + A(AX + iAY) = UZ + \AZ. (2.19) 

Consider the evolution equation for the process Z^(t, x) = X^{t, x) + zy'^(t, x), 

di 

dZ\t,x) = Y [Lr,xZ\t, ■) + Fr{t,x)] dV; + Y [SiZ\t,x) + Gi{t,x)]dMl, 

r=0 1>1 

(2.20) 

Z\0, x) = Zo{x) = Xo{x) + iYo{x). (2.21) 

In order to prove well-posedeness of the problem (I2.20l) - fl2.21l) . firstly we have to 
check the following stochastic parabolicity condition: 

Condition (CI) There exists a constant K >0 such that for Z e t e [0,T]: 

di 

2j2{LrZ , Z)m dV[ + Y \\Si{Z) \\l d{M')t < K\\Z\\l dVt 

r=0 «>0 

a.s. in the sense of measures. 

Let Z = X + e using ^Tm and fl^TTD . we deduce that 

di 5 

2 Y ' ^"^) + 5Z 5Z ^(^')* = ^^-W' (2.22) 

|a|=m r=0 \a\=m 1>1 k=1 

where to ease notation we drop the time index in the right handsides and we set: 

d di 



\a\=mj,k=l r=l 

+ {al'^'D^D'^Y, DkD'^Y) + {V^^DjD'^X, DkD'^Y) 
+ J2 \{(^jDjD''X,afDkD''X) + {ajDjD''Y,a^DkD''Y) d{M^) 



dV[ 



Z>1 



dT2{t) 



di 
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r=0 |a|<m j,k=l (/3,7)e2:(a):|7|=l 



+ 



^{D^ai'''DjD^X,DkD'^X) + {D^ai'^DjD^Y, DkD'^Y) 

d 

[{aiDjD''X,D''X) + {aiDjD''Y,D''Y) 



dT3{t)=Y,Y.H E [{D^'criD,D^X,afD,D-X) 

1>1 j,k=l \a\=m (/3,7)G2(a);|7| = l 

+ {D^'ajD,D"Y,afDkD''Y)d{M')t, 
d 

dn{t)=J2Yl E [{^iDjD''X,aiD-X)-iajD^D"X,riD»Y) 

j,k=l \a\=m 

+ {aiDjD''Y,TiD''X) + a^D^r)] d(M')t, 



\a\V\l3\<m r=0 j,k=l 
d 



+ J2{C'X)D^Z,D"Z) + {Cr{.)D^Z,D''Z) 



dv: 



a 

+ E [E ( [C'KO + E Ci{-)]]D^Z, D-Z)] d(M'),}, 
«>1 j=l 1>1 

where C^'*^, C^, C^, C/ and Q are bounded functions from M°' to C due to Assumption 
(A4(m,p)) for any p E [2, oo). 



For every r the matrix bj. is symmetric; hence ^, [(6^' i^ji^^X, D^D^'Y) — 
{bi.'^DjD''Y,DkD''X)] = 0. Hence, Assumption (A2) used with yj = DjO'^X 
and with yj = DjD°'Y, j = I,-- - ,d, imphes Ti{t) < for t G [0,T]. Further- 
more, Assumption (A3(m)) yields the existence of a constant C > such that 
dT^it) < C\\Z{t)\\'^dVt for all t e [0,T]. Integration by parts shows that for regular 
enough functions f,g,h:W^^ R, 7) G X(a) with I7I = 1, we have 



ifD^g,D"h) = -{fD"g,D^h) - {DyD^g.D^h). 



(2.23) 



Therefore, the symmetry of the matrices implies that for (j) G {X{t)^Y{t)} and 
r = 0, ■ ■ ■ 



j,k=i 



d 



j,k=i 
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A similar argument proves that for fixed j = 1, ■ ■ ■ ,d, r = 0, ■ ■ ■ ,di and = X(t) 
or = Y{t), 

{aiDjD''(f),D"4>) = -^{DjaiD"(P,D''(f)). 

Therefore, Assumption (A3(0)) imphes the existence of if > such that dT2{t) < 
K\\Z{t)\\'^dVt for all t G [0,T]. Furthermore, dT/^it) is the sum of terms {(pip, Djcj)) 
d{M^)t where = D"X{t) or cp = D''Y{t), and ip = fg, with / E {af} and 
g G {ai,Ti}. The identity {<pip,Dj(p) = —^{(pDjip,(p), which is easily deduced from 
integration by parts, and Assumptions (Al) and (A3(m)) imply the existence of 
K >0 such that dTi{t) < K\\Z{t)\\l^dVt for every t G [0,T]. The term dT^it) is the 
sum over I > 1 and multi-indices a with |a| = m of 

d 

j,k=l (/3,7)eX(a):|7| = l 

with = X{t) OT = Y{t) and // = aj for every j = 1, ■ ■ ■ ,d. Then, A[l, a) = 
B{1, a) — C{1, a), where B{1, a) = Yl'j k=i ^j,k{^^ Q^) and 

(/3,7)e2(a):|7|=l 
d 

C{1, a) = J2 E [DV^fiD.DP^, D.D'^^) . 

i,fc=l(/3,7)gX(Q);|7|=l 

Integrating by parts twice and exchanging the partial derivatives Dj and Dk in each 
term of the sum in C(/, a), we deduce that 

{Dytf^D,D^V,D,D\) = -{[D,[DytfnDjD^V + DytfiDkD,D^^],D'^^) 
= ( - Du[Dy^fi]D,DPip + {DAD^ftmD.Df^^, 

On the other hand, by symmetry we obviously have 

j,k j,k 

Using Assumptions (Al) and (A3(l)) we deduce that there exist bounded functions 
(p{a, a, I) defined for multi-indices a which have at most one component different 
from those of a, and such that 

A{l,a) = ^B{l,a) + ^ a, OD'^'I', . 

\a\=m 
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Furthermore, integration by parts yields 

d 1 
j,k=l j,fc=l (/3,7)GX(a):|7l=l 

Thus, Assumption (A3(l)) imphes the existence of a constant C > such that for 
the various choices of ^ and /f, Ei>iE\a\=m\B{l,(y)\d{M')t < C\\Z{t)\\ldVt for 
every t G [0,T]. Therefore, we deduce that we can find a constant A' > such that 
dT^{t) < K\\Z{t)\\'^dVt. The above inequahties and fl2.22p complete the proof of 
Condition (CI). 

Since Lr are linear operators. Condition (CI) implies the following classical Mono- 
tonicity, Coercivity and Hemicontinuity: for every Z,( ^ jjrn+i defined by 

din, 

di 

2Y,{LrZ -LrC,Z- OmdV; + ^ \\Si{Z) - Si{0\\ld{M% < K\\Z - C\\ldVt, 

r=0 1>1 

2j2{Lr,xZ , Z)^dV; + E \\Si{Z)\\ld{M), + 2X\\Z\\l^,J2v: < K\\Z\\ldV, 

r=0 1>1 r=0 

a.s. in the sense of measures, and for Zi G H"^~^^, i = 1,2,3, r = 0, ■ ■ ■ ,di and 
A > 0, the map a G C — )■ {Lr^x{Zi + aZ2 , Zs)m is continuous. 

The following condition (C2) gathers some useful bounds on the operators Lr 
and Si for < r < di and / > 1. 

Condition (C2) There exist positive constants Ki,i = 2,3,4 such that for Z G 
A G [0, 1], r = 0, ■ ■ ■ , di and / > 1: 

2\\Lr,xZ\\l^, + \\SiZ\\l < K2\\Z\\l^, a.s. (2.24) 
\iSiZ,Z)m\ < Ksll^ll^and \{SiZ,Gi)m\ < K4Z\\rn\\Gi\U+i a.s.. (2.25) 

The inequality (12.241) is a straightforward consequence of the Cauchy-Schwarz in- 
equality and of Assumption (A3(m)). Using integration by parts and Assumptions 
(A3(m))-(A4(m,p)), we deduce that if Gi{t) = Gz,i(t) + iG;,2(i), 

d 

\{SiZ,ZU <- E + {D,alD'^Y,D'^Y)]\ 

\a\=m j=l 

+ Yl \{C{a,P,l)D^Z,D^Z)\ 

a|V|/3|<m 

d 

\{SiZ, GiU < E E (W^"^' DjD^Gi^r) I + I [alD^Y, D.D^Gi^^) \ 

j=l \a\=m 
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+ Yl \{C{a,P,l)D^Z,D^Gi))\, 

\a\V\[S\<m 

for constants C{a, (3, 1), a, /3, 1 such that sup^,^^ ; C(a, f3, 1) < C < oo. Hence a simple 
apphcation of the Cauchy-Schwarz and Young inequahties imphes inequahty f l2.25p . 

We then proceed as in the proof of Theorem 3.1 in [TB] for fixed A > (see also 
j2Uj and [ID])- To ease notations, we do not write the Galerkin approximation as 
the following estimates would be valid with constants which do not depend on the 
dimension of the Galerkin approximation, and hence would still be true for the weak 
and weak* limit in L'^{[0,T] x Q; Hm+i) and L'^^Q; L°°(0,T; Hm)). Let us fix a real 
number iV > and let tn = infjt > : ||^^(t)iU > N} AT. The Ito formula, the 
stochastic parabolicity condition (CI) and the Davies inequality imply that for any 
t e [0,T] and A G (0,1], 



E sup \\Z\s A T^)\\l) + 2AE / \\Z\s)\\l^,ds < E\\Zo\t 

^ s&[0,t] ^ Jo 

+ 2Ve/ \{Fr{sATN).Z\sATr,))m\dV: 

+ 5^E / [2\{SiZ\s ATr,),Gi{s AT^))J + \\Gi{s AT^)\\l]d{M'), 

l>l 

+ 6 e( ^ { r {SiZ\s A Tn) + Gi{s A rjv), Z\s A r^))^(i(M'),}' 

1>1 

The Cauchy-Schwarz inequality, the upper estimate ( ]2.25p in Condition (C2) and 
inequalities (12. 7p - (12. 8p in Assumption (Al) imply the existence of some constant 
A' > such that for any 5 > 0, 

E sup ||Z^(sAr^)||^, +2AE / \\ZXs)\\l^,ds < E\\Zo\t 
^ se[o,t] ^ Jo 

+ 3 5e( sup \\Z\sATM.)+K5-^y,^ / \\Fr{s)\\ldV: 
^se[o,t] ' 7^ Jo 

t rt 

\t „ A _ M|2 



+ e/ 5^K[ri + l]||GKs)||^^+irf(M'). + E / \\z> 
Jo Jo 



[sATM)\\'dVs. 



For (5 = |, the Gronwall Lemma implies that for some constant G we have for all 
AT > and A G (0,1], 

E sup ||Z^(.Ar^)||^ +AE / \\Z\s)\\l^,ds < GE{\\Z4l + KUT)) . 
^se\o,i\ ^ Jo 
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As — oo, we deduce that r^r — t- oo a.s. and by the monotone convergence 
theorem, 

cT 



Ef sup \\Z\s)\\l) +AE 
^ se[o,T] ^ 



\Z^(s) 



_,ds<CE{\\Zo\\l + KUT)). 



Furthermore, Z'^ belongs a.s. to C([0, T], H"^). As in [TB], we deduce the existence of 
a sequence A„, — ?■ such that Z^'^ — )■ Z weakly in L^([0,T] x ^2; H"^). Furthermore, 
Z is a solution to fl2.ip and f l2.2p such that (12.1 5p holds and is a.s. weakly continuous 
from [0, T] to H"^. 

The uniqueness of the solution follows from the growth condition fl2.24p in (C2) 
and the monotonicity condition which is deduced from the stochastic parabolicity 
property (CI). 

(ii) Suppose that Assumption (A4(m,p)) holds for p G [2, oo). Set p = 2p 
with p G [l,oo); the Ito formula, the stochastic parabolicity condition (CI), the 
growth conditions (C2), the Burkholder-Davies-Gundy and Schwarz inequalities 
yield the existence of some constant Cp which also depends in the various constants 
in Assumptions (Al)-(A4(m,p)), and conditions (C1)-(C2), such that: 



E( sup 

se[o,t] 



<cJe||Zo||^ + e(sup \\Z{s)\\i 



E 

r-=0 



\FJs)\ir.dv: 



E sup \\Z{s)rm 



E sup \\Z[s) 



I m 



E 

i>i 

E 



\Gi{s)U^^d{M'), 



Zis) 



+ \\Gi{s)\\l]d{M^), 



p/2 



Using the Holder and Young inequalities, (I2.13P as well as Assumptions (Al) we 
deduce the existence of a constant A' > such that for any 5 > 



e( sup \\Z{s)\\l) < 35e( sup \\Z{s)\\l\+C{p)\. 



z[0 

+ AfP^-^E 



t "1 



Eii'-^ 



idv: 



\nz,rm 



r=0 

t 



1>1 



Let = I and introduce the stopping time rjy = inf{t > : > A^} A T. 

Replacing t hj t A in the above upper estimates, the Gronwall Lemma and 



fl2.13p prove that there exists a constant C such that E^ sup^gjQ j]^^^ ||Z(s)||^j < 

CE(||Zo||^ + A'.m(T)^/^) for every A^ > 0. As A^ — t- oo the monotone convergence 
theorem concludes the proof of fl2.16p . This ends of the proof of Theorem [1] □ 
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Remark 2.1. If ai'^{t,x) = 0, for example for the Schrddinger equation, Assump- 
tion (A2) implies that aj = 0. 

2.2. Further a priori estimates on the difference. Theorem [1] is used to upper 
estimate moments of the difference of two processes solutions to equations of type 
(lOl . Fore = 0,1, r = 0,--- , di,l>l,j,k = !,■■■ ,d,\et ai'^^{t,x),bl^^^{t),ai^{t,x), 
ae^r(t,x),b^^r(t,x),ali(t,x),as,i(t,x),Ts,i(t,x) be coefficients, F^^r(t,x), Ge,i(t,x) be 
processes, and let Z^q be random variables which satisfy the assumptions (Al)- 
(A3(m)) and (A4(m,p)) for some m > 1, p G [2, oo), the same martingales 
{Ml,t G [0,T]) and increasing processes {y[ ,t G [0,T]). Let L^^^ and S^^i be de- 
fined as in (12. 3p and (12. 4p respectively. Extend these above coefficients, operators, 
processes and random variables to e G [0, 1] as follows: if /o and /i are given, for 
e G [0, 1], let fe = efi + (1 — £)fo- Note that by convexity, all the previous assump- 
tions are satisfied for any e G [0, 1]. Given e G [0, 1], let denote the solution to 
the evolution equation: Zs{0,x) = Zi^^o^x) and 

di 

dZ,(t,x) = YiLs,rZeit,x)+F,^rit,x)]dV[+Y, [Ss,iZsit, x)+Gs,iit, x)]dMl (2.26) 

r=0 l>l 

Thus, Theorem [T] immediatly yields the following 

Corollary 1. With the notations above, the solution Z^ to (I2.26P with the initial con- 
dition Zi^fi exists and is unique with trajectories m C{[0, T]; H"'"^) f] L~(0, T; H"^). 
Furthermore, the trajectories of Z^ belong a.s. to Ctu([0, T]; if™) and there exists a 
constant Cp > such that 

sup e( sup ||^.(t,-)ll^) <Cp sup ¥.[\\Z,^^f^^ K^{TYl''\ < oo. (2.27) 

eG[0,l] ^tG[0,T] ^ £G{0,1} ^ ' 

Following the arguments in [TU], this enables us to estimate moments of Z^ — Zq 
in terms of a process which is a formal derivative of Z^ with respect to e. Given 
operators or processes f^, s G {0, 1}, set /' = /i — /o- 

Theorem 2. Let m > 3, and p G [2, oo). Then for any integer k = 0, - ■ ■ , m — 2 

e( sup ||Zi(t)-Zo(t)||^) < sup e( sup WUtWX (2.28) 

where is the unique solution to the following linear evolution equation: 

di 

dC{t) = Y,{Ls,rCeit,x) + L'^Z,{t,x) + F;.{t,x))dVr{t) 
r=0 

+ J2 {SsAit, x) + S'iZ,{t, x) + G'lit, x))dMl , (2.29) 
i>i 

with the initial condition Z'^ = Zi — Zq. Furthermore, 

sup e( sup ||C.(t)||^_2) <oo- (2.30) 
eG[o,i] ^te[o,r] ^ 
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Proof. Using f l2.27p we deduce that the processes Fr(t,x) = Lj.Z^(t, x) + F^(t, x) and 
Si{t,x) = SiZi;{t,x) + G'i{t,x) satisfy the assumption (A4(m-2,p)) with m — 2 > 1. 
Hence the existence and uniqueness of the process Ce, solution to (12.291) . as well as 
(I2.30p can be deduced from Theorem [1] 

We now prove (12.281) for k, G {0, ■ ■ ■ , m — 2} and assume that the right hand-side 
is finite. Given (/e, e G [0, 1)), for and h > and e E [0, 1] such that e + h E [0, 1], 
set 6hfe = ife+h — fe)/h. We at first prove that (12.281) can be deduced from the 
following: for every e G [0, 1), as /i — )■ is such that h + e E [0, 1], 



E 



( sup \\5hZe{t)-Ut)\\i) ^0. (2.31) 



Indeed, assume that (I2.3ip holds and for n > let i?„ = n'^A^(nld — A)"'' denote 
the fcapa-fold composition of the resolvent of the Laplace operator A on the space 
l? = . Then, by some classical estimates, there exists a constant C{k) > such 
that for any (p G L"^, WRnhW^ < C(k)||</)||o. Hence (12.311) yields that for every n > 0, 

as /i with e + h e [0, 1], we have E^sup^gjo.r] \\5hRnZe{t) - RnCe{t)V^^ -> 0. 

Furthermore, since for every integer > 1, we have Zi — Zq = J2k=o ^i/NZk/N < 
supjg[o 1] 6i/]\fZir, we deduce that for every n > and p G [2, oo): 

e( sup \\RnZo{t) - RnZ,{t)\\l) < sup e( sup \\RnCe 

Finally, if G is such that hminf^^oo ||-Rn0lU — < oo, then (p G H'^ and 
II 011 K < Thus, by applying the Fatou Lemma and using estimate (12.301) we can 
conclude the proof of (I2.28p . 

We will now prove the convergence (I2.3ip . It is easy to see that the process 
Ve,h(t, ■) := 5hZe{t, ■) — Ce(^, ") has initial condition ?7e,/i(0) = 0, and is a solution of 
the evolution equation: 



r=0 

+ Yl [Se,lVe,h{t, ■) + Sl{Z,+f,{t, ■) - Z,{t, ■))]dM^ 



1>1 



Hence, using once more Theorem [T], we deduce the existence of a constant Cp > 
independent of e G [0, 1) and h > 0, such that e + /i G [0, 1], 



sup WShZeit) - Cemi) < CpE( [ ||Z,+,(t) - Z,it)\\ldVt 
efo.Tl ^ ^Jo 



p/2 



1 /3 2 /3 

Using the interpolation inequality ||0||2 < C'H^Hq H^Hs , see for instance Propo- 
sition 2.3 in [17], the Holder inequality and the estimate (12.270 with m = 3 from 
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Corollary [H we deduce that 



e( sup \\ShZ,{t) - C{t)\\l) <c\e( sup \\Z,+hit) - Z,{t)\\l 

Vt6[0,T] ^ L \/crnTl 



te[o,T] 



1/3 



Finally, the process ^e,h{t, ■) = Z^+hit, ■) — Zs{t, ■) is solution to the evolution equa- 
tion 

di 
r=0 

+ Yl [Se,Mt, ■) + hS'iZ.^hit, ■) + hG[{t, ■))]dMl, 
l>l 

with the initial condition ^^ /^(O) = h{Zi — Zq). Thus, fl2.27l) and Theorem [T] prove 
the existence of a constant C, which does not depend on 5 G [0, 1) and h > Q with 
e + h E [0, 1], and such that 

Ef sup \\Z,+,,{t) - Z,{t)\\i) <CK''\ 
^te[o,r] ^ 

This concludes the proof of (12.311) and hence that of the Theorem [2l □ □ 

3. Speed of convergence 

3.1. Convergence for time-independent coefficients. For r = 0, ■ ■ ■ , c^i, e = 
0,1, let (yj^,t G [0,T]) be increasing processes which satisfy Assumptions (Al), 
(A2), (A3(m-|-3)) and (A4(m-|-3,p)) for some integer m > 1, some p G [2, +00) 
separately for the increasing processes (y^^,t G [0,T]), the same increasing process 
{Vt,t G [0,T]) and the initial conditions Z^ q, e = 0,1. For e = 0,1, let Z^ denote 
the solution to the evolution equation 

dZ,{t,x) = [LrZe{t,x) + Fr.{x)] dV^, + ^ [SiZe{t,x) + G^(x)]c^M,^ (3.1) 

0<r<di l>l 

with the initial conditions Zq{0, ■) = Zq q and Zi{0, ■) = Zi^ respectively. Let 

A := sup sup max — V'q jI- 

i^en t€[0,T] r=0,l,- -,di 

Then the H'^ norm of the difference Zi — Zq can be estimated in terms of A as 
follows when the coefficients of and Fr are time-independent. Indeed, unlike the 
statements in [12] , but as it is clear from the proof, the diffusion coefficients ai and 
Gi can depend on time. 

Theorem 3. Let Lr and Fr be time-independent, J^Q-measurable, VJ, e = 0,1, Mi 
be as above and let Assumptions (Al), (A2), (A3(m-|-3)) and (A4(m-|-3,p)) be 

satisfied for some m > and some p G [2, +00). Suppose furthermore that 

E(|X;i|i^r|tL+/ + sup \Y\\Gris)\\l_,/ )<00. (3.2) 



r=0 
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Then there exists a constant C > 0, which only depends on d and the constants in the 
above assumptions, such that the solutions Zq and Zi to (13.1 p satisfy the following 
inequality: 

e( sup ||Zi(t)-Zo(t)||^) <c(e(||Zi,o-Zo,o|D + Ap). 

The proof of Theorem [3] will require several steps. Some of them do not depend on 
the fact that the coefficients are time independent; we are keeping general coefficients 
whenever this is possible. The first step is to use Theorem [2] and hence to define a 
process Z^ for any e G [0, 1]; it does not depend on the fact that the coefficients are 
time-independent and extends to the setting of the previous section. For e G [0, 1], 
r = 0, ■ ■ ■ rfi, t G [0, T] and x G M^, let 

Vl, = eVl, + il-e)V,:,, pl, = dVlJdVt 

and for j,k = set ai'^{t,x) = pl fii'^{t,x), bi'^{t) = pl,.bl'''{t), ai^^{t,x) = 

Pl,tK{t^^)^ a£,r(t,a;) = pli^arit.x), K,r{t,x) = pljbrit.x), L^^r = Plt^r, F^,r{t, x) = 

plj.Fr{t,x). Then for e G [0,1], the solution Zsit,-) to equation (12.11) with the 
increasing processes VJ'^ can be rewritten as (I2.26P with the initial data ^e(O) = 
eZifi + (1 — e)ZQfi and the operators (resp. processes) S^a = Si (resp. Ge^i = Gi). 
Furthermore, we have 

di d 

Y.^'^'i^e!;it,x)+^bi'\t))dv: = 

r=0 j,k=l 

d 

0<r<di j,k=l 

Hence the conditions (Al), (A2), (A3(m-|-3)) and (A4(m-|-3,p)) are satisfied. 
Therefore, using Theorem [21 one deduces that the proof of Theorem [3] reduces to 
check that 

sup e( sup ||C.(t)||^) <C{E\\Z^^o-Zo,oK+AP), (3.3) 
£e[o,i] ^te[o,i] ^ 

where if one lets = — V^^, the process Ce is the unique solution to (I2.29P 
which here can be written as follows: for t G [0, T] one has 

di 

dCeit) = J2 [^'■^-(^' ^^"^^It + {LrZeit, x) + F,(t, x))dAl] 
r=0 

+ J2SiCe{t,x)dMi{t), (3.4) 

i>i 

and the initial condition is Ce(0) = ^i,o ~ ^o.o- 

To ease notations, given a multi-index a, j,k E {!,■ ■ ■ ,d} and Z smooth enough, 
set Z^ = D'^Z, Zo,,j = D'^DjZ and Zo,,j,k = D'^D^D/.Z, so that for Z, C G H"", 
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{Z, C)m = Y.\a<m(Za, Cq)o- Let 

^ = I ^^a^'^Z^Z^ ; a""'^ uniformly bounded and complex- valued, Z G 

a p 

and for $, ^ G ^ set $ ~ if there exists Z G H"^ such that - '^){x)dx = 

Jjgd r{x)dx, where F is a function defined by 

r(x) = J2 Zc,{x)P^Z{x) with = ^ 7"'^^/3, (3.5) 

for some complex- valued functions 7"'^ such that |7"'^| are estimated from above by 
the constants appearing in Assumptions (Al), (A2), (A3(m-|-3)), (A4(m-|-3,p)). 

Note that if F is as above, then for some constant Cm(r) we have 

/ \T{x)\dx<C^{T)\\Z\\l. (3.6) 
For e > 0, J, A; = 1, ■ ■ ■ d, / > 1 and t G [0, T], set = d{M^)t/dVt and let 



r=0 ^ «>1 

For m > and z G H"^^^, set 

[^]™ := = E («f W^.^, D,Z)^ + Cr^\\Z\\l, (3.7) 

with Co = and Cm. > to be chosen later on, so that the right handside of (13. 7p 
is non negative. Given Z, C G H^+\ set [Z, C]m = \{[Z + Ql + [Z - C]V)- We at 
first prove that defines a non-negative quadratic form on for some large 

enough constant Cm- Once more, this result does not require that the coefficients 
be time-independent. 

Lemma 1. Suppose that the conditions (Al), (A2) and (A3(m-|-1)) are satisfied. 
Then there exists a large enough constant Cm such that (13. 7p defines a non-negative 
quadratic form on . 

Proof. Assumption (A2) for e G {0, 1} implies that (13. 7p holds for m = and 
Co = 0. Let m > 1 and a be a multi-index such that 1 < |a| < m. The Leibnitz 
formula implies the existence of constants C(q;,/3,7) such that for Z G H"^~^^, 
d 

5^(af(t)D,Z)^,Z„,,)^ = 

j,k=i 

d 

J2 (sf (t)^«, , ^a,fc)o + E c{a,p,^)i:'^'^{t), 

j,k=l /3+7=a,l/3|>l 
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where I^'^''^{t) := Yl'^jk=i ' -^a,fc)o' Furthermore given m > 1, multi- 

indices a with |a| < m and Z G H"^~^^, using Assumption (A2) we deduce that 
Si<jfc<d {H'^i^)^a,j 5 ^a,k)Q > 0. Thus, the proof reduces to check that 

I^'^'\t) ~ 0. (3.8) 

Indeed, then the upper estimate (13.61) proves (13. 7p . which concludes the proof of the 
Lemma. Integration by parts implies /"■^''^(t) = — Xl^fc=i {j^k[^D^ay^{t)Z^ j) , . 
Since |/3| < m, by Assumption (A3(m-|-1)) we know that DkD^a{'^{t) is bounded 
and hence /,"'^'^(t) ~ - Elfc=i(^''«P(^)^7,i,fc > ^«)o- l/^l ^ 2, then I7I < m - 2; 
hence by (A3(m)) we deduce that /^"''^•^(t) ~ 0. If = 1, then ap(t) = fi^^(t), 
so that 

I ^''5f(t,x)D^(X,,,(.)X,,,(.) + K,,(.)K,,,(.))(x)rfx. 

Thus, integrating by parts and using (A3(2)) and the inequality I7I + 1 < m, we 

deduce that 4"'^'^(t) ^J2lk=i {D^D^H'Hi)Z'y,j ^ ^7,fc)o ~ 0- ^his concludes 

the proof. □ 

The following lemma gathers some technical results which again hold for time- 
dependent coefficients. 

Lemma 2. Suppose that the assumptions of TheoremlMhold. There exists a constant 
C such that for ( G H^~^^ one has dVt a.e. 

p(c):=2 J2 pm(c,^.c)^ + E^*ii^'^ii™ + 2[C<^^iiciil (3.9) 

0<r<di 1>1 

For any f = 0,1, ■■■ ,di, Z e H"^+^ and ( G let 

QfiC, Z)= ^ Pl,t [{LrC, LfZ).m + (C, LfLrZ)ra] + ^ qliSlC, LfSlZ)m. 
0<r<di />1 

Then there exists a constant C such that for any Z G and C, G if ™"*"^, one has 

dVt a.e. 

|g,(C,^)l<C||Z|U+3(||CIU,+ [C]m). (3.10) 

Proof. Suppose at first that C ^ i/™""*"^; since the upper estimates (13. 9 p and (I3.10p 
only involve the if'""'~^-norm of C,, they will follow by approximation. Then we have 

2p;,(C,L,C)™ + 5^gl||5,C||^,= $^g?(c,c), 

0<r<di l>l \a\<m 
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where 

gr(C,C) = 2 Pit {CaALrOo.), + Y.<it\\{SiOX- 

0<r<di 1>1 

Integration by parts and assumption (A3(m)) imply that for |a| < m, one has 
2(Ca , KrQ)a)o ~ 2 / a^^it, x) {X^{x)X^j{x) + Y^{x)Y^,,{x))dx 
= / al/t,x){Xl + Y^)(x)dxr^-{al,it),Ca,Ca),-0, 

- - I (a/aO),(x)|G(x)|2rfx~0, 
||(K + zrfc)C)«||'~0. 

Finally, we have (Ca , (A(^&^.''=(t)/^,C),)^ = 0. Set L% = E^',^, (apD,C) 
and SiC, = Yl'j=i i'^i + ^'^/)-^jC- Then we have 

Qr(C,C)~2X;PM(C«, (^°C)a)o + E^*ll(^°Oallo- (3-11) 

r=0 Z>1 

If m = 0, integration by parts proves that the right hand side of (13. lip is equal to 
— 2[(^]q (with Co = 0). Let m > 1 and a be a multi index such that m > \a\ > 1; 
set r(a) = {(/3,7) : a = /3 + 7, |/3| = 1}. For 0, e H"", let C(/3,7) be coefficients 
such that: 

D''{(j)ij) = (PD^'^/j + C{(3,-f)D^(f)D^tlj + Y C{(3,-f)D^(f)D^^lj. 

(/3,7)er(a) /3+7=a,|/3|>2 

This yields 

E "'J II WO jiM: "! E { K w ■ c-)o 

1>1 1>1 j,k=l 

(/3,7)er(a) 
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Since for (/3,7) G T{a) we have I7I + 1 < |a| < m while + 1 = 2, integrating by 
parts and using (A3(m)) we have for fixed /, 

j,k j,k 

Furthermore, integration by parts and (A3(m)) yield 

- E C^(/3,7)(^^K^)C7,,,C«,fc)o}- 

(/3,7)Gr(a) 

Therefore, the definition of a^'j:, ( I3.1ip and ( 13. 8 p yield 

hk j,k (/3,7)er(a) 

--2E((o«f).^c...)^. 

Hence for C G 

^'(c) = E / ^"(^' + ml = 2 E ^"c)o, (3.12) 

a|<m. |o|<m 

for some operator which satisfies f l3.5p . Hence fl3.6p concludes the proof of fl3.9p . 
Polarizing f l3.12p . we deduce that for ZX ^ H"^^^, 

di 

J2 Pit [{Z, LrOm + {LrZ, Qm] +Y,^t {SlZ , SlC)m. + 2[Z , Qm 
r=0 1>1 

\a\<in 

Let f e {0, 1, ■ ■ ■ , rfi} and for Z e C e set Z = LfZ; then if one sets 

di 

QfiC^Z) := Yp'^^t[{LrC,LfZ)m + {C,LfLrZ)m] +E^* {SiCLfSiZ)^, 

r=0 i>l 

one deduces that 

q,iC,z)+J2p:,t{C^[LrU-UU]z)^+2j2ql{SiC [SiU-L,Si]z)^+2[c,uzu 

r I 

= E [{D"UZ,P"Oo + {D%,P^UZ)o\. 

\a\<m 
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The operators L^Lf — LfL^ and SiLf — LfSi are of order 3 and 2 respectively. Hence 
integration by parts and the Cauchy Schwarz inequahty imply that 

\qr{C,Z)\<C\\an\\Z\\rn+3 + C[CUZU. 

Finally, (13. 7p and Assumption (A3(m)) imply that for Z G H"^~^^, 

< C'll^llm+i + C'mll^llm < C'll^llm+1- 

This concludes the proof of (13.101) . □ 

The following lemma is based on some time integration by parts and requires the 
coefficients of Lj. and Fj. to be time independent. 

Lemma 3. Let the assumptions of Theorem\^ he satisfied and Zg, (resp. Q) denote 
the processes defined by (13. ip (resp. (13. 4p ). For r = 0, ■ ■ ■ ,di and t G [0,T], let 
K = Vlt - V^, and set 

Js,t := J2 f (^-(^)' ^^^^(^) + ^r)jA:. (3.13) 
Then there exists a constant C such that for any stopping time r <T , 

E sup (Je,t- / [Us)fmdVs) 
l-te[0,r] ^ Jo ^ + ^ 

< sup \\Us)\\1)+c(a^ + ¥. fwUsWrndv). (3.14) 

Proof. The main problem is to upper estimate J^^t in terms of A and not in terms 
of the total variation of the measures dA\. This requires some integration by parts; 
equations (13 .ip and (13. 4p imply: 

di 

J,, = (Ut) , UZ,{t) + Fr)^Al - J2 Je,t, (3.15) 

r=0 l<fc<4 

where for t G [0,T] we set: 



Jit =J2^^s [ X^(^rC.(s), LrZ,{s) + Fr)m + (C.(s) , (s) + F,]), 

r f 

Jl = r E E {Siis)Us). LASi{s)Z,{s) + Gi{s)])j{M^),, 

-^0 r 1>1 

+ (C.(s),L,[^K^)^.(s) + G'Ks)])Jt/Mi, 



dAl. 
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Note that 

r,f 

Using f l3.10p . integration by parts, Assumption (A3(m)), the Cauchy-Schwarz and 
Young inequahties, we deduce that 



Jit + Jit <CA \ \\ZM\U,{[US)U + \\Ce{s)\\m} + J2 11^^(^)11- 

+ CA [ y2\\Ceis)\\^\\Giis)\Usd{M')s 
Jo I 

< f\[Ce{s)]l+\\Cs{s)\\l)dV, 



+ CA^ 



\ZJs 



sK-^J llm+3 



< {\US)]1 + WUsWm) dVs + CA'^I^ 11^.(^)11^+3 + Km+2{t) 

where the last inequahty is deduced from Assumption (A4(m+2,2)). 

The Cauchy Schwarz inequahty, integration by parts and Assumption (A3(m+1)) 

imply that for fixed r = 0, ■ ■ ■ ,di and / > 1, 

I {S^s)Us), LrZ^is) + F,.)J + I (C.(.), Lr[Siis)Z,is) + Giis)])J 

< C\\Ce{s)\\m[\\Ze{s)\\m+3 + ||Fr|U+l + ||GKs)|U+2]. 



Therefore, the Burkholder Davies Gundy inequality and Assumption (Al) imply 
that for any stopping time r < T, we have: 

e( sup |Jir/2 



m+3 ' m+l 



0<r<di 

p/4 



sup Ell^'(^)ll-+2 rf(M'). 

S6[0,T] J 



< CAP/^E 



di 



sup \\Zs{s)ti + \j2m 

X {l^\\cs{s)\\idv:)"^'] 



2 

m+l 



p/4 



sup 5^11^^(5)11^+2 

SG[0,T] I 



p/4 



1 
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sup 



se[o,T] 



m+3 



+ \Kj2\\Pr\\ 



2 

m+1 



r=0 



p/2 



sup VllGK^) 

selo,T] ' ^ 



lm+2 



«>1 



p/2' 



+ -E sup iic.(.)ii!;)+cE / imwi^dVs. 

op ^se[o,T] 



Using the condition (13.21) and Theorem [T] with m + 3, we deduce yhat 

e( sup IJJJP/^) < ^e( sup \\Ut)rm)+CE r \\C{s)\\ldVs + CAP. 



Therefore, 



e{ sup Ut- [\Ce{s)Us)dV,Y'] < 1/(8p)e( sup IICWII 

+CE r \\Ce{s)r^dVs+cAP+cE( sup 

Jo ^ t&[0,T] 

Integrating by parts we obtain 



2Ji = J2 {LfUt) + F, , LrZ^it) + Fr)^A\JV, - J] jj/, (3.16) 



where 



r,f 1>1 •^'^ 

Jtf = E E r ^'sK {L,[S,{s)ZM + Gi{s)] , L,[SK.)^.(s) + Gz(5)])^rf(M'),. 
„^ ,^1 ~'o 



Integration by part. Assumption (A3(m+2)), the Cauchy-Schwarz and Young in- 
equalities yield 



[Lr[LrZ,{s) + Fr] , L,Z,(s) + F,) J < C[||^.(s) 



I m+3 m+2 
|2 



\{Lr\S,{s)Z,{s) + G,{s)],LrZ,{s) + F,)J<C[\\ZM\\^^^^ 
Hence, using Theorem [1], (13.21) . Assumptions (Al) and (A4(m+2)) we deduce 

E( sup \j^:l + J^f\'^')<CA^. 

Finally, the Burkholder-Davies-Gundy inequality implies that 
e( sup 14^/^ 
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< CA^E 



{Lr[Siis)Z,{s) + Gi{s)] , LrZ,{s) + Fr)Jd{M'), 



p/4 



< CA^. 



Hence, E(supjg[o,r] \ Je,tf^'^) — CA^, which concludes the proof. 



□ 



Using Lemmas [T]13l we now prove Theorem [3] for time-independent coefficients. 
Proof of Theorem Apply the operator D°' to both sides of fl3.4l) and use the Ito 
formula for ||D"(^£(t)||Q. This yields 

|a|<m r 
|a|<m 1>1 

where {Z, Qm denotes the duality between and which extends the scalar 

product in H^. Using (13. 9p we deduce that 

dWUml < -2Ut)]ldVt + CWUmldVt + 2dJ,, + 2^(Ce(t), Si{t)Ut))jMl, 

1>1 

where Je,t is defined by fIXTSD . Using (12:251) we deduce that \ {Ceit), Si{t)C{t)) J < 
C\\(e{t)\\m- Thus Lemma [H the Burkholder Davies Gundy inequality and Assump- 
tion (Al) yield for any stopping time t <T 



e( sup WCeitWm) <CE\\Zr,o-Zo,orm+P^( SUp Je,t - [ [Ce{s)]ldV, 



+c,E / wcemtdVs 



p/4 



+ CpE 



[0/ 

\\Ceis)\\ldVs 



p/2 



p/2 



< 



CE||Zi,o-^o,o||L + ^e( sup \\Ut)rm)+c(A^ + E r WUsWmdVs) 
4 ^te[o,T] ^ ^ Jo ^ 

C,E rWUsWmdVs 

Jo 

ce( r\\us)r^dv. 



+ CE 



sup WCeitWrlH I \Us)\tdV. 

te[o,r] 



p/4- 



< CE||Zi,o - ZoflVm + CAP + Ie( sup lie 

^ ^t£[0,T] 



p 

'£ llm 



where the last upper estimate follows from the Young inequality. 

Let Tat = inf{t : ||Ce(^)|lm — ^} ^ ^! then the Gronwall Lemma implies that 

e( sup ||C.(t)L) <C(E||Zi,o-^,o||i; + ^^). 
Letting — > oo concludes the proof. 



□ 



ON THE SPLITTING METHOD 25 

3.2. Case of the time dependent coefficients. In this section, we prove a con- 
vergence result similar to that in Theorem [3] when the coefficients of the operators 
depend on time. Integration by parts in Lemma El will give extra terms, which 
require more assumptions to be dealt with. 

Assumption (A5(m)) There exists an integer number d2, an (j-'t) -continuous 
martingale A''^ = {N^, • ■ ■ , N^^) and, for each 7 = 0, ■ • • , 0^2 a bounded predicable 
process : Q x (0, T] x M'^ — )■ M.^ for some N depending on d and di such that 

h^{t, X) := (a^;t(^' 3;), «7,r(^' «7,r(^, X), b^,r{t, x) , x)] 

l<j,k<d, <r <di, 1 < 7 < rfs), 

for some symmetric non negative matrices {aii^^{t,x),j,k = 1, ■ ■ ■ ,d) and 
j,k = 1, - ■ ■ ,d). Furthermore, we suppose that for every u E Q and t G [0, T], the 
maps h^{t, ■) are of class C""*"^ such that for some constant K we have \D°'h^{t, ■)\ < 
K for any multi-index a with |a| < m + 1 and such that for t G [0, T], 

7=1 

h{t,x) = h{0,x)+ / ho{s,x)dV, + y2 h^{s,x)dN]. 

Jo 7=1 >^o 

For 7 = 0, ■ ■ ■ , (^2, T = 0, ■ ■ ■ , (ii, let be the time dependent differential operator 
defined by: 

d d 

L^,rZ{t,x) = A-(K',t(^'^) + iV^Ut)]DjZ{t,x)^ + J2<rit,x)D,Z{t,x) 

j,k=l j=l 

+ [a^^rit, x) + ib^^rit, x)~\ Z(t, x). 
For r = 0, ■ ■ ■ ,di, let 

LrZ{t, x) = Lr{0)Z{0, x)+ Lo^rZ{s, x)dVs + y] L^,rZ{s, x)dN], 

and Fr{t,x) = Fr{0,x) + Fo,^(s, + E7L1 ^7,r-(s, a;) c^iVJ. We then have the 
following abstract convergence result which extends TheorerrO 

Theorem 4. Let Assumptions (A(l)), (A(2)), (A3(m-|-3)), (A4(m-|-3,p)) and 
(A5(m)) be satisfied and suppose that 

e( sup I J2 + sup I 5^ l|G„+2(t)||^,+2 < 00. (3.17) 

^t€[0,T]' tcrn-rll'' — 



r=0 



te[o,T] I 
such th( 

e( sup \\Z,{t)-Zo{t)rj<c(E{\\Z,{0)-Zo{OWrn)+AA. (3.18) 



Then there exists some constant C > such that 
sup 

•t&[0,T] 
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Proof. Since Lemmas [T] and [2] did not depend on the fact that the coefficients are 
time-independent, only Lemma [3] has to be extended. For t G [0,T], let 

Since = for r = 0, ■ ■ ■ ,di, the integration by parts formula (13.151) has to be 
replaced by 

di 7 
Je,t = i^eit) , Lr{t)Z,{t) + Fr{t))^Al - Je,U t G [0, T], 

r=0 k=l 

where the additional terms on the right hand-side are defined for t G [0, T] as follows: 

jl,=Y,A:{Ceis),Lo,rZM + Fo,r)jVs, 
r 

Jl = fY.A:f2{Cs{s),L,,rZe{s) + K,As))jNl 
Jo ^ 

"^0 r 1>1 7=1 



7= 

Arguments similar to those used in the proof of Lemma [3l using integration by parts 
and the regularity assumptions of the coefficients, prove that for k = 5,6 there exists 
a constant C > such that for any stopping time tau < T we have: 

t€[0,T] ' ^te[0,T] tG[0,- 



e( sup \Jl,\^l^\<CA^I''^[ sup WmVl^ sup (||Z,(t)|U^2 + C)' 
<-^e( sup \\at)€)^CA\ 



where the last inequality follows from the Young inequality. Furthermore, the 
Burkholder-Davies-Gundy inequality and the upper estimates of the quadratic vari- 
ations of the martingales A^"^ yield for every 7 = I,-- - ,^2, r = O,--- ,di and 
r G [0,T], 

sup |JJ,r/2^ r \A:{Us),L,,rZs{s) + F,As))JdVsY' 

I'r!' sup {\\ZMU2 + CY^"' 

te[o,T] ' 
Hence, the proof will completed by extending the upper estimate fl3.16p as follows: 



E 



<CA^/2E sup Wamt^ sup (||Z,(t)|U^2 + C)' 

^tg[0,r] i6[0,r] 



7 
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where for j = 4, ■ ■ ■ , 7 we have: 

r,f "^0 

- ^0 

r,r 7,7 

Jef = 2 E E E r ^^'^^ (^7,K^)^.(^) + F,As) , 
„^ ,^1 Jo 



r,f 7 ;>1 



Je'J = 2 E E /* (^7,r-(5)^.(^) + F,,,{S) , Lr{s)Z,{s) + Fr{s)) JN]. 

^, Jo 



r,r 7 

T4,fc| 



We obtain upper estimates of the terms E( sup^^jQ^] | J^'^ 1^/^) for /c = 4, ■ ■ ■ , 7 by ar- 
guments similar to that used for /c = 1, ■ ■ ■ ,3, which implies E( sup^^jQ l^e^J^^^) < 
CA^. This concludes the proof. □ 

4. Speed of convergence for the spliting method 

The aim of this section is to show how the abstract convergence results obtained in 
Section [3] yield the convergence of a splitting method and extends the corresponding 
results from [TU]. The proof, which is very similar to that in [TU] is briefly sketched 
for the reader's convenience. 

Assumption (A) Forr = 0, ■ ■ ■ , di, let L^- he defined by (12.31) and for I > 1 let Si be 
defined by fl2.4l) . Suppose that the Assumptions (A2) and (A3(m-|-3)) are satisfied, 
and that for every u G Q andr, I, Fr(t) = Fr(t, ■) is a weakly continuous H"^~^^-valued 
function and Gi{t) = Gi(t, ■) is a weakly continuous H'^'^'^ -valued function. Suppose 
furthermore that the initial condition Zq G L^{Q;H"^'^^) is J^Q-measurable, that Fr 
and Gi are predictable and that for some constant K one has 

di 

e( sup Ell^^(^'-)ll-+3+ sup J]||G(t,-)||L+4+ll^oL+3) 

Let = (y^^,t G [0,T]) be a predictable, continuous increasing process such that 
Vq = and that there exists a constant K such that + ^;>i(Af')r < K . Finally 
suppose that the following stochastic parabolicity condition holds. 
For every {t, x) G [0, T] x M.'^, every uj ^ fl and every A G M^, 



E 

j,k=l " 1>1 



> 



in the sense of measures on [0,T]. 
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Let Z be the process solution to the evolution equation: 

dZ{t,x) =(LoZ(t,x) + Fo(t,x))rfrj° + ^ {LrZit.x) ^ Fr{t,x))dt 

r=\ 

\^{SiZ{t,x)\Gi{t,x))dM\ (4.1) 

with the initial condition 2^(0,-) = Zq. Then Theorem [T] proves the existence and 
uniquness of the solution to (14.11) . and that 

e( sup II Z(t) 11^+3) <C 

for some constant C which depends only on rf, di, m, p and T. 

For every integer > 1 let 7^ = {ti := iT/n,i = 0, 1,-- - ,?7,} denote a grid 
on the interval [0,T] with constant mesh S = T/n. For > 1, let Z^"'^ denote 
the approximation of Z defined for t G 7^ using the following splitting method: 
Z(")n(0) = and for z = 0, 1, ■ ■ ■ , n - 1, let 

Z^-\U^^) := ) ■ ■ ■ Pf )g,^,^,,Z(")(t,), (4.2) 

where for r = 1, ■ ■ ■ ,di and t G [0, T], P/'^V denotes the solution (r of the evolution 
equation 

d(r(t,x) = (^LrCr(t,x) + Fr(t,x))dt and (^r(0,x) = ipi^x), 

and for s G [0, t] < T, Qs,t4' denotes the solution 77 of the evolution equation defined 
on [s,T] by the "initial" condition ri{s,x) = ip{x) and for t G [s,T] by: 

dr]{t,x) = {Lor]{t,x) + Fo{t,x))dVt^ + ^{Si7]{t,x) + Gi{t,x))dMl. 

i>i 

The following theorem gives the speed of convergence of this approximation. 

Theorem 5. Let a^''', b^'^, al, a,., br, crj, Fr, Gi satisfy the Assumption (A) . Sup- 
pose that al'^,bl'^,ai,ar,br and Fr are time-independent. Then there exists a con- 
stant C > such that 

E(^^ ||Z(")(t) - Z(t)ll^) < Cn~P, for every n > 1. 
teTn 

Proof. Let d' = di + 1 and let us introduce the following time change: 

f for t < 0, 

K{t) = { t- k6di for t G [kd'6, {kd' + 1)5), A; = 0, 1, ■ ■ ■ , - 1, 

y {k + 1)5 for t G [{kd' + 1)5, {k + l)d'5), A; = 0, 1, ■ ■ ■ , n - 1. 

Let, for evry t G [0, T], 

it) = M^(,) , = , = V,', = y° , 

y^y^ = K(t), v;y^ = ^(t - r(5) for r = 1, 2, ■ ■ ■ , di. 
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For e = 0, 1, consider the evolution equations with the same initial condition 
Zq{0,x) = Zi{0,x) = Zq{x) and 

di 

dZ,{t) = {LrZe{t) + Fr)d% + + Gi)dMl. (4.3) 

r=0 

One easily checks that the Assumptions (Al), (A2), (A3(m+3)), (A4(m+3,p)) 

are satisfied with the martingales Mi and the increasing processes V^j for 5 = 0, 1 
and r = 0, 1, ■ ■ ■ , c?!. Therefore, Theorem [1] implies that for e = 0, 1, the equation 
(14.31) has a unique solution. Furthermore, since condition (13.21) holds. Theorem [3] 
proves the existence of a constant C such that 

Ef sup \\Ziit) - ZQ{t)\\A < C sup max I fi:(t + r5) -K(t)|P = CTPn"P. 

Since by construction, we have Z^^d't) = Z{t) and Zi{d't) = Z^'^\t) for t G 7^, this 
concludes the proof. □ 

Note that the above theorem yields a splitting method for the following linear 
Schrodinger equation on W^: 

d 

dZ{t, x) = (^iAZ{t, x) + Y tt' {x)DjZ{t, x) + F{x) ] dt 

+ Yl [(^'(^) + ^nix))Z{t,x) + Gi{x)]dMl, 
i>i 

where a\ F (resp. ai, ti and Gi) belong to H^^^ (resp. H"^~^^). Indeed, this model 
is obtained with a^''^ = aj = and IP''' = 1 for j,k = 1, - ■ ■ ,d and I > 1. 

Finally, Theorem H] yields the following theorem for the splitting method in the 
case of time demendent coefficients. The proof, similar to that of Theorem O will 
be omitted; see also [TU], Theorem 5.2. 

Theorem 6. Let ai''', bl''', ai, aj, ai, ti, F, Gi satisfy Assumptions (A) and (A5(m)) . 

For every integer n > 1 let Z*^"^ be defined by (14.21) when the operators L^, Si, 
the processes F^ and Gi depend on time in a predictable way. Then there exists a 
constant C > such that for every n>l, we have: 

E(5^||Z("Ht)-^(t)r™) <Cn-^ 
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